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Abstract Iron based narrow gap semiconductors such as FeSi, FeSb2, or FeGa3 
have received a lot of attention because they exhibit a large thermopower, as well 
as striking similarities to heavy fermion Kondo insulators. Many proposals have 
been advanced, however, lacking quantitative methodologies applied to this prob- 
lem, a consensus remained elusive to date. Here, we employ realistic many-body 
calculations to elucidate the impact of electronic correlation effects on FeSi. Our 
methodology accounts for all substantial anomalies observed in FeSi: the metal- 
lization, the lack of conservation of spectral weight in optical spectroscopy, and 
the Curie susceptibility. In particular we find a very good agreement for the anoma- 
lous thermoelectric power. Validated by this congruence with experiment, we further 
discuss a new physical picture of the microscopic nature of the insulator-to-metal 
crossover. Indeed, we find the suppression of the Seebeck coefficient to be driven by 
correlation induced incoherence. Finally, we compare FeSi to its iso-structural and 
iso-electronic homologue RuSi, and predict that partially substituted Fei_ A Ru v Si 
will exhibit an increased thermopower at intermediate temperatures. 
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1 Introduction 

Correlated semiconductors have been a subject of intensive research over the 
years, because they exhibit an unusual metalization process which is poorly un- 
derstood. At low temperatures, the iron silicide FeSi - the prototypical compound 
of this class of materials - is akin to an ordinary semiconductor with a gap of 
A « 50 — 60meV JT] |2] [3] |4j- Yet, at higher temperatures, that are however much 
smaller than A /kg, FeSi becomes a bad metal|5j |2] |4| and develops a Curie-Weiss 
like susceptibility HI |6). Analogies with heavy fermion Kondo insulators [8] [9] 
and mixed valence compounds IflOl , effects of spin-fluctuations lfTTI . or spin-state 
transitions lfl2l [T3ll , as well as an anomalous electron-phonon coupling lfPfl [TBI [T6l 
have been invoked to account for this behaviour. Lacking, however, quantitative 
methodologies applied to this problem, a consensus remained elusive to date. 

Besides this fundamental puzzle, the class of correlated narrow gap semiconduc- 
tors, which also comprises compounds such as FeSb2, FeAs2 or FeGa3, is also of 
interest in view of applied science. Indeed these materials exhibit notably large See- 
beck coefficients at low temperatures 131 [T71 [181 [TBI [191 |20l |2T1 |22l . and the largest 
thermoelectric powerfactor ever to be measured was recently found for FeSb? lfT8l . 
An understanding of the electronic structure of these compounds, and the potential 
influence of electronic correlation effects thereon is thus of vital interest. 

Here, we will investigate the paradigmatic example FeSi by means of the realis- 
tic extension of dynamical mean-field theory (see e.g. Ref. ||23l for a review). Our 
recent results l24l account for all substantial anomalies observed in FeSi, namely the 
lack of conservation of spectral weight in the optical conductivity l2l [3l l25l[T3l l4l l26ll . 
a Curie-Weiss like susceptibility and an anomalous thermoelectric power. From the 
microscopic insight of our approach, we elucidate the origin of the metal-insulator 
transition. We explain that the latter is a consequence of a correlation induced in- 
coherence: Unlike in conventional semiconductors where a metalization process is 
driven by thermal activation or a moving of the chemical potential into the conduc- 
tion or valence bands, in correlated insulators such as FeSi the crossover is induced 
by the emergence of non-quasiparticle incoherent states in the gap. We further link 
the occurrence of these many-body states to the spin degrees of freedom. 

With this understanding we here address the Seebeck coefficient of FeSi, and 
propose ways on how to improve the thermoelectric performance of FeSi based 
systems. In particular, we compare FeSi to its isoelectronic homologue RuSi, and 
speculate on the thermoelectric properties of Ru substituted FeSi. 
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Table 1 Band-structure of FeSi in the B20 structure for different degrees of the distortion from the 
simple cubic rock-salt structure. Energies in electronvolts, the origin corresponds to the Fermi level. 
From left to right: "NaCl" rock-salt structure 0%, 50%, 75%, and 100% distortion corresponding 
to actual FeSi. Above the band-structures are shown the respective unit cells (pictures made with 
1271 ). iron atoms are red, silicon ones blue. 



2 The metal-insulator transition in FeSi and its microscopic 
origin 



2.1 Crystal- and band-structure 



The B20 crystal structure of FeSi, despite being cubic (space group P2i3), is rather 
complex and has four iron atoms per unit cell. However, it can be viewed as a sim- 
ple cubic rock-salt structure that is highly distorted along the [111] direction[280- 
In order to elucidate the important effect of this distortion onto the band-structure, 
we display in Fig. Q~]the evolution of the Kohn-Sham spectrum (using the gener- 
alized gradient approximation (GGA) functional within density functional theory 
(DFT) as implemented in wien2k |29l ) for different degrees of the atomic displace- 
mentsj. While FeSi is metallic in the fictitious rock-salt structure, the changes in 
the atomic positions cause the formation of an avoided crossing in the distorted 
B20 structure. In Fesi, the iron 3d orbitals split into a low-lying z 2 , and the two 
doubly degenerate groups x 2 -y 2 , xy and xz, yz, which are separated by a band 
gap lF28l 1X41 l30l l24l . With a nominal valence of 6 electrons per iron, band-structure 



1 Indeed, the positions of the iron and silicon atoms are (u,u,u), (j+u, j-u,l-u), (l-u,^+u,j+u), 
and ( |-u,l-u, j+u) with u(Fe)=l/4 and u(Si)=3/4 for the rock-salt structure, and u(Fe)=0.136 and 
u(Sn=0.844 for FeSil28l. 

2 as measured by a linear interpolation of the u parameters from above (see also Ref. 1 30 1). 
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methods thus correctly reproduce the low temperature insulating character of FeSi. 
Owing to the hybridization nature of the gap, as well as the absence of many-body 
band-narrowing effects within DFT methods, the gap is overestimated by a factor of 
two, A DFT «0.11eVES[H|30|. 
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Fig. 1 Spectral function. LDA+DMFT spectrum at various temperatures. The inset shows a zoom 
of the low-energy pseudo-gap, and grey symbols indicate the spectral minimum, which moves as 
a function of temperature owing to the particle/hole asymmetry of the spectrum. 



2.2 The many -body spectrum 

While the shown band-structure may qualitatively describe the low temperature 
properties of FeSi, e.g. the activated behaviour found in the resistivity below 
150K|5][17), it obviously cannot explain the insulator to metal transition at higher 
temperatures, let alone a Curie-Weiss like tail in the spin susceptibility. 

To account for electronic correlation effects beyond the effective one-particle 
description of band-theory, we resort to the realistic extension "LDA+DMFT" of 
dynamical mean-field theory l23l . as implemented in Ref. 1311 . We use established 
values of the local Coulomb interaction in iron based compounds (32 (U=5.0eV, 
J=0.7eV), and solve the DMFT impurity by means of a continuous time quantum 
Monte Carlo (ctqmc) method I331 I 351 . 

In Fig. [T] are displayed the local spectral functions that we obtain for different 
temperatures as indicated. The low temperature spectrum is akin to the density of 
states obtained within band-theory, yet with a gap that is renormalized by about a 
factor of two, in agreement with photoemission spectroscopy experiments 11351 [36ll . 
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This means that lifetime effects are minor at low temperatures, as also inferred from 
the self-energy shown in Fig. [2] In particular the spectral weight at the edges of the 
gap is remarkably spiky, a characteristic that is commonly considered a hallmark for 
a potentially large thermopower ll37l . With increasing temperature however, spectral 
features broaden, and the charge gap gets filled with incoherent weight, with solely 
a pseudo-gap remaining. This is in congruence with experimental findings, both in 
one-particle probes (such as photoemission[35 36 1), transport measurements (for a 
comparison of our theoretical resistivity with experiments, see Ref. 124)), and opti- 
cal spectroscopy (again, see Ref. J24)). Besides the filling of the gap, a further detail 
of photoemission results is captured: Since the spectrum is particle/hole asymmet- 
ric, the chemical potential moves as a function of temperature. As a measure for this, 
we mark in the inset of Fig.Q]the position in energy at which spectral weight is min- 
imal. At low temperatures, the latter is in the vicinity of the chemical potential (as 
expected for a semiconductor|21|). Since the unoccupied states have larger spectral 
weight than the valance states, the chemical potential moves down upon increasing 
temperature, causing the spectral minimum to move up in energy, as indeed found 
in photoemission l35l l36l . Above 300K, the asymmetry switches, and the position 
of the spectral minimum moves back towards the chemical potential. These trends 
in the particle/hole symmetry are harbingers of the temperature dependence of the 
Seebeck coefficient, discussed in Section|3] 




50 100 150 200 250 300 350 400 
T[K] 

Fig. 2 Correlation induced incoherence. Imaginary part of the DMFT self-energy at the Fermi 
level as a function of temperature and resolved into orbital characters. Also shown is the quadratic 
fit: 3X(o = 0) = - 1.9 ■ 10~ 4 meV/K 2 ■ T 2 for the x 2 -y 2 lxy and xzlyz components. 

Thus, as far as spectral properties are concerned our calculations are a realistic 
generalization of the seminal model of Fu and Doniach l38l . However, as we detail 
in Ref. 11241 . there is a fundamental physical difference between that model and our 
results. While the degree of correlations in the two band model f38l is controlled by 
Hubbard physics, we find our results to be much more sensitive to the Hund's rule 
coupling J than the Hubbard U. This highlights the multiorbital nature of the system 
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and, hence, the necessity of realistic calculations. For a more detailed discussion, 
see Ref. ll24l . In the context of FeSi note also the recent Ref.j39), as well as Refs. 
I40ll4l"ll42l for the influence of the Hund's rule coupling in other systems. 



2.3 The physical picture - the self-energy and the spin state 

As seen above, the insulator to (bad) metal transition is not caused by a closure of 
the charge gap, but by a filling of the latter with incoherent spectral weight. The 
degree of incoherence of the many-body system is encoded in the imaginary part of 
the electron self-energy, SE(co). In Fig. [2] we depict the zero frequency limit of that 
quantity, resolved into the iron 3d orbital characters. The orbital components that 
account for spectral weight at low energies are, as mentioned above, the x 2 -y 2 and 
xy on the valence, and the xz and yz on the conduction side. As seen in the figure, the 
latter two neatly follow a T 2 law for the temperature range shown, with a coefficient 
of - 1 .9 • l(T 4 meV/K 2 . Therewith the inverse lifetime quickly becomes comparable 
to the size of the gap, which it surpasses at about 400K. 

Having thus ascribed the metalization process in FeSi to an effect of electronic 
correlations, we may ask as to the physical origin of this coherence-decoherence 
crossover. For this we look at the microscopic insights furnished by our theoretical 
approach. Within dynamical mean field theory the system is described by an ef- 
fective impurity that represents the iron atoms, and a hybridization function (Weiss 
field) that accounts for the embedding of the reference system into the solid. Useful 
information can be gained by decomposing the local projection of the system into 
the eigenstates of the effective impurity, as shown in Fig. [3] Displayed is the proba- 
bility of the local reference system to be in a state with N particles and spin state 5, 
for (a) low, and (b) high temperature. 

The histogram shows a rather broad distribution over many spin and charge 
states, signaling large fluctuations at short time scales. Indeed the variance of the 
charge state 8N = ((N — (N)) 2 ) ~ 0.93 is of order unity, entailing an overall mixed 
valence state, with an average 3d occupation (N) = 6.2. Also virtual spin fluctua- 
tions are large, 55 = ((5 — (5)) ) « 0.33. Further, we obtain from the spin distribu- 
tion an effective moment M = y/S(S+ l)g s ~ 3 (g s = 2), which is consistent with 
major contributions from 5=1 states, and in congruence with the moment M = 2.7 
as obtained from fitting a Curie-Weiss law % = ^j^-M 2 / (T — T c ) to the experimen- 
tal susceptibility fT| [6] for T > 400K. The dominance of states with 5=1 implies 
in particular that FeSi is not a singlet insulator lTT2l [131 . 

Further, we note that the decomposition shown in Fig. [3] (a) & (b) is basically 
independent of temperature. This means that the mixed valence state is not tem- 
perature induced, as previously proposed lfTol . Moreover, this also rules out a spin 
state transition. This is to be contrasted to systems like MnO, or LaCo03 in which 
high-spin/low-spin transitions occur l43l l44l l45l l46l . 

The temperature independence of the variances of the spin and charge state 
(short time fluctuations) is in stark contrast to the experimental uniform suscepti- 
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(b) 
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Fig. 3 The spin and charge state. Decomposition of the DMFT impurity onto the eigenstates of 
the effective atom, resolved into the spin S and charge state N. 



bility (time averaged response) which shows activated behaviour at low T, and a 
Curie-Weiss like decay at temperatures beyond 400K. Thus, there is a tremendous 
differentiation in time-scales: While the underlying spin structure of the (effective) 
iron atoms does not evolve, the manifestation of the fluctuating moment in the spin 
response is highly susceptible to external conditional 

The temperature induced unlocking of the fluctuating moment of the iron sites 
establishes a link to real space. Therewith the momentum space description of the 
low temperature coherent semiconductor breaks down and effective lifetimes are 
introduced as the system decomposes over states of different momenta. 



3 The Seebeck coefficient of FeSi and RuSi 

Finally we turn to the discussion of the thermoelectric properties of FeSi. For details 
of the employed linear response formalism, see Refs. [21 , 24], and e.g. Refs l47ll48l 
49 50]|5T|[52| for related prior works. 

Before addressing the influence of many-body effects onto the thermopower, we 
note that it was found that -for temperatures below 100K- the Seebeck coefficient of 
FeSi can actually be reproduced by a slightly hole doped band-structure [ 14|. This 
we confirm by adding 0.001 holes per iron as well as an overall effective mass of two 
to the density functional results from above. The thus obtained Seebeck coefficient 
indeed yields a good agreement at low temperatures, as shown in Fig. |4] In our 
opinion, the hole doping should not be viewed as an introduction of extra charge, 
but as a way to adjust the particle/hole asymmetry of the band-structural 



3 For the theoretical local spin susceptibility see Ref . 1 24 1 . 

4 However, a notable dependence of the thermopower on the precise stoichiometry is witnessed in 
experiments ll7l . 
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200 400 600 800 1000 1200 
T[K] 

Fig. 4 The Seebeck coefficient of FeSi and RuSi. While a description of FeSi requires employing 
a many-body theory, the thermopower of RuSi is reproducible based on band-structure methods, 
see text for details. Experimental data from [5 15 1 (FeSi) and |53| (RuSi). Indicated with thin 
dotted lines are the maximal electronic contributions to the Seebeck coefficient A /T for both com- 
pounds. 

The congruence of a band-structure based thermopower with experiment fur- 
ther corroborates that FeSi is a coherent semiconductor at low temperatures. In 
the regime kgT <C A the thermopower of such a system is - modulo the temper- 
ature independent Heikes contribution - given by (A/2 ~ jx)/T ■ 8X, where A is the 
charge gap, fi is the chemical potential measured from the centre of the gap, and 8X 
quantifies the electron/hole asymmetry (for details see Ref. |2T1 ). Since the latter 
is constrained, \8X\ < 1 {8X = +1 would e.g. correspond to a purely hole driven 
thermopower), the electronic contributions to the Seebeck coefficient are, for an in- 
sulator, limited to ±A /T. Thus, assuming the same asymmetry, a larger gap causes 
a greater Seebeck coefficient. As indicated in Fig.|H FeSi respects the above bound- 
ary, while it is, for example, largely surpassed in the related compound FeSb? |[2T|| . 
advocating for that material the importance of non-electronic contributions, espe- 
cially the phonon-drag mechanism. The quenching of the Seebeck coefficient at 
very low temperatures (S — > for T — > 0) is the consequence of a small yet finite 
scattering rate ETlPi 

The use of a renormalized band-structure fails at describing the properties of 
FeSi above 100K, when the metalization process becomes notable, as seen in the 
above Fig. Q] for the spectral function (photoemission experiments), and now in 
Fig.|4]for the Seebeck coefficient. Extending our scheme 12T1 for the thermopower 

5 For transport calculations that are based on density functional theory we assume a constant scat- 
tering rate 3X = — lmeV. 



Thermopower of the Correlated Narrow Gap Semiconductor FeSi 



9 



to include dynamical self-energy effects, we calculate the Seebeck coefficient based 
on our realistic dynamical mean field results. As displayed in Fig. [4] the thus ob- 
tained thermopower is in very good agreement with experimental findings. In par- 
ticular, we capture the changes of sign as a function of temperature which indicate 
the transition between hole (S > 0) and electron (S < 0) dominated transport. The 
non-monotonous tendencies in the Seebeck coefficient were already heralded by 
the moving of the chemical potential, as seen in the inset of Fig. [T] Starting from 
low temperatures, the chemical potential moves down, therewith reducing the hole 
contributions to the Seebeck coefficient, before it passes, at around 120K, the point 
of thermoelectric particle/hole symmetry, below which the Seebeck coefficient be- 
comes negative. At yet higher temperatures, the trend reverses and the thermopower 
changes sign again. 

Thus, in the current case correlation effects are detrimental for the thermoelectric 
performance. This has to be contrasted to the case of correlated metals, where a 
reduced quasi-particle weight increases the Seebeck coefficient (see e.g. Refs. l48l 

E2). 

We find it instructive to compare FeSi to its iso-structural and iso-electronic ho- 
mologue RuSi, and discuss the series Fei_ r Ru A Si in view of its potential as thermo- 
electricfl Ruthenium silicide, RuSi, is a semiconductor with a gap of 0.2 — 0.3eV 
as inferred from optical spectroscopy l55l l56l or resistivity measurements ll55l l53l . 
Interestingly the tendency in the size of the gap in Fei_ A Ru v Si is not monotonous 
in xj57). Indeed up to a ruthenium concentration of 6% the charge gap is found 
to decrease. With increasing concentration the gap then grows, bypasses the initial 
value of FeSi at around 15% ruthenium, and further augments up to the value of 
pure RuSi. 

With our insight into FeSi, we can explain this trend. As a matter of fact, there 
are two opposing tendencies that have to be considered: On the one hand, ruthe- 
nium has a larger atomic radius than iron. Thus with increasing Ru content, the 
hybridization gap (see Section 12. Q will shrink as the lattice expands. This effect 
is immediate, and wins for low ruthenium concentrations. On the other hand, one 
has to consider the crossover in the dominant orbital character of excitations at low 
energy, namely the transition from 3d to 4d electrons. In FeSi we identified the ratio 
of the Hund's rule coupling J and the bandwidth as the controlling factor for the 
strength of correlations l24l . Since the 4d electrons of ruthenium are less localized 
than the iron 3d ones, this ratio decreases with growing ruthenium admixtures, re- 
sulting in smaller effective masses and thus a weaker many-body narrowing of the 
charge gap. This effect becomes preponderant for larger ruthenium concentrations. 

This decrease in electronic correlation effects, of course, also means that the 
use of conventional band-structure methods is more justified for RuSi than for FeSi: 
Using again the GGA functional within wien2k l29l . we obtain a band-structure (not 
shown) that resembles that of FeSi, albeit with stronger dispersions and a larger gap. 
The latter is, as in FeSi, indirect, and amounts to 0.23eV, in agreement with previous 
band-structure calculations [[58, 59, .60 ] . For the Seebeck coefficient, we here find for 



6 For a comparison of the related couple FeSb2 and RuSbi, see Refs. 1201 1541 . 
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RuSi - as was the case for FeSi - that a small hole doping is needed to adjust the 
particle/hole asymmetry in density functional results. Indeed, as shown in Fig. 21 
doping pure RuSi with only 0.0025 holes per ruthenium, yields a thermopower in 
excellent agreement with experiment l53l . 

The Seebeck coefficient of RuSi, while not reaching the very large values of 
FeSi at low temperatures, is notable in size over an extended temperature regime, 
with about 250/xV/K from 100 to 500K. Comparing RuSi with FeSi, we make two 
interesting observations: (a) as noted above, the thermopower of a coherent insulator 
is controlled by the size of the gap and the particle/hole asymmetry. RuSi having a 
larger gap than FeSi, its Seebeck coefficient indeed surpasses the envelope function 
of FeSi, A FeSi /T, for 200K < T < 500K, as indicated in Fig. g] However, it is 
further away from its own boundary, Ar^/T, than is the case for FeSi. Hence, 
the particle/hole asymmetry (above called 8X) is smaller in the 4d compound, (b) 
Besides the peak value which is smaller in RuSi, the temperature dependence is 
qualitatively akin. Indeed, when scaling the temperature with the respective sizes of 
their low temperature gaps, the curves for FeSi and RuSi are very similaiQ. 

Without actually performing many-body calculations for RuSi, we can expect 
that also in other experimental observables the temperature dependence scales with 
the ratio of the respective gap values. In particular the emergence of a fluctuating 
moment that causes the coherence-decoherence crossover and quenches the Seebeck 
coefficient in FeSi will be pushed to higher temperatures. 

Therefore, we believe that the iso-valent substitution of iron with ruthenium, 
Fei_ x Ru x Si, will provide competitive thermoelectric properties in the temperature 
regime of 100 to 250K. Besides the greater coherence, the loss of thermoelectric 
particle/hole asymmetry (8X) when departing from pure FeSi, is partly compensated 
by the enlarged gap A (for x > 0.15). While the power factor S 2 G, with the dc 
conductivity a, of pure FeSi reaches favourable 40 /iW/(K 2 cm) at around 60K l24l . 
we expect a notable improvement over both pure FeSi and pure RuSi beyond 100K. 
Finally, the substitution will also reduce the thermal conductivity k, yielding a better 
figure of merit ZT = S 2 (JT/k. 



Conclusions 

In conclusion, we presented a new scenario for the intriguing properties of iron 
silicide FeSi, in which the metalization process with temperature is driven by corre- 
lation induced incoherence that we traced back to the unlocking of fluctuating iron 
moments. Using realistic many-body techniques, we investigated the signatures of 
this microscopic theory in the one-particle spectrum, the optical conductivity (see 
Ref. ll24l ). and the Seebeck coefficient, for all of which we find quantitative agree- 
ment with experiment. 



7 It would be of great value to have experimental measurements on RuSi up to higher temperatures. 
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With the thus gained physical insight into the interplay of electronic correlation 
effects and thermoelectricity, the general class of narrow gap semiconductors be- 
comes more amendable to theory assisted thermoelectric material design. Indeed, 
we made explicit suggestions for improvements of FeSi based thermoelectrics. 
While pure FeSi already presents a competitive power factor in the temperature 
range of 50-100K, we proposed means for extending the favourably large Seebeck 
coefficient to higher temperatures, by motivating a study of Fei-jRu^Si vis-a-vis its 
thermoelectric performance. 
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